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Abstract 

The magnetized relativistic Fermi and Bose gases are studied at finite tem- 



lO perature and density. In the case of the Fermi gas, the contribution to the 

f — . magnetization from the vacuum becomes dominant for high magnetic fields, 



when the thermal contribution saturates. In the case of the charged Bose gas, 

^ the (paramagnetic) vacuum-magnetization becomes dominant when the gas 

changes from a diamagnetic to a paramagnetic behaviour. We furthermore 

find that the scalar-QFD effective coupling constant for a weak non-zero ex- 



<D ternal magnetic field is a decreasing function of the temperature. 
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1 Introduction 

The magnetostatic properties of the vacuum has hitherto not been studied to the same 
extent as the electrostatic case where the vacuum polarization has well-known screening 
effects. However, it is necessary to add the vacuum contribution to the magnetization in 
order that the average of the induced current equals the curl of the magnetization (see 
Ref. [f]). Then also the mean field equation follows from the minimization of the effective 
action with respect to the magnetic induction B (see below). The vacuum magnetization 
becomes dominant in spinor- as well as scalar-QED, at extremely high field strengths. 
The magnetic fields associated with collapsed magnetic stellar objects may be as high as 
B = 0(fO^°)T [2] (cf. mj/e = O{10^)T). At such field strengths the vacuum magneti- 
zation starts competing with the ordinary thermal magnetization. Here we study some 
aspects of a an ideal QED e"'"e~-plasma and a charged Bose gas at finite temperature 
T and chemical potential fi in the presence of a static uniform magnetic field, in partic- 
ularly the high field limit. The thermal Fermi and Bose gases in a high magnetic field 
have been studied earlier in Refs.[3, 4], but there the vacuum contribution was neglected. 
Since our numerical algorithms are converging faster at vanishing chemical potential, we 
shall mostly consider the neutral plasma, but the general behaviour is unchanged at finite 
density. 

2 The Fermi Gas 

The partition function Zj(i?,r, //) for the relativistic ideal e"'"e~-plasma in presence of 
an external magnetic field i? in a sufficiently large quantization volume V, and the corre- 
sponding thermal part of the effective Lagrangian jC^JI^ , can be written as [5] 
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Here the energy spectrum is -Ea,™ = y^^ + 2eB{n -|- A — f ) -|- m^, and fp{E) are the 
Fermi-Dirac equilibrium distributions 



exp[/3(^ T /«)] + 1 



where /3 is the inverse temperature. Separating the field independent part, we write 

£f/ = £^'^ + £f^, where 

° " m7 " 3^ i-oo ^^^ ~ ^ )/^(^) (^ - "1 j . (3) 

Here Zf[T^fi) is the partition function for the field-free ideal e"'"e~-plasma with particle 
energy E = \/k^ + m^ and fpi'^) = ^(^)/f (^) + ^(~^)/f (~^)- Using the identity 

exp( — |a;|) [^ dt f ^ ^ 2 , ^ 
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jC^,'" = ^J2i-^y^l — exp -^^ m^x cosh(/3//i)[e5xcoth(e5x)-l] .(5) 



exp --(x^t + -) , (4) 



/o V27rt V 2' t 

and expanding the distribution functions under the restriction |/i| < m, we can write the 
field dependent part Cl'^ for |/i| < ?7i as 

In the case // = 0, Eq.(5) agrees with the result obtained in Ref. [6] where it was derived 
using thermal quantum field theory. Here we have shown that it can be directly derived 
from the canonical partition function. It is, however, not obvious how to generalize Cl'^ 
to |/i| > ?7i, since then it appears to be divergent due to the presence of the infinite sum in 
Eq.(5). The latter problem has recently been solved and we refer to Ref. [5] for a general 
discussion. The effective Lagrangian density is 

£efr = /:o+/:i+/:f/ , (6) 



where the tree-level part is Cq = — I/2i?^, and C\ corresponds to the the well-known 
result [7] 

[eBy [^ dx / 2 ^ \ f 1 2I 

£1 = ^-r — exp(-m ^)<^a;coth(a;) - I - -x ^ . (7) 
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We have here performed the standard renormalizations leaving eB invariant. Notice that 
Eq.(7) corresponds to a renormalized / = term in Eq.(5). We then obtain 
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where we have identified a ??4-function, given by ??4[2;, (j'] = 1 + 2 X]/^i(~l)'<?' cos(2/2;). The 
effective Lagrangian can be used to define an effective temperature dependent coupling 
constant [8] 
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For weak fields, eB <^ m^ , and for T ^ m tlie last term in Eq.(8) dominates and the 
corresponding effective coupling agrees with a conventional renormalization group calcu- 
lation [5], which can be showed using the expansion J2'hLi Ko{xl){ — l) "'"^ — )• —1/2 log x as 
X ^0. 

An external field is included by adding a term C^xt = jext ■ A to CeS- Here jext is the 
external current which is independent of the dynamics of the system considered, such 
that V X H = jext, and A is the vector potential B = V X A . Neglecting a boundary 
term at infinity we rewrite C^xt = B • H. The mean field equation 

B = H + M{B) , (iO) 

is then obtained by minimizing the effective action with respect to i?, where we have 
included the vacuum contribution Afyac in the magnetization, i.e. 

M = M'^'^ + M,,, = ^(£f'^ + £i) . (if) 

We find that the thermal part M^'^ of the magnetization saturates for high fields (a = 

eV47r) 



eM^'" ^ — T^ ,V7B>T->m . (i2) 

This was discussed in Ref.[3], but without considering the vacuum contribution 

eMvac ?^ — e51og(e5/m2) , eB :> m^ , (i3) 

that starts to dominate when eBlog[eB /m^) > tt^T^. We have numerically evaluated 
Afvac and M^'^ and the result is presented in Fig.f. In order to improve the convergence 
in the numerical calculation we have found it convenient to write Eq.(8) in the form 
J^ dxf[x) = JQdxf{x) + J^ dxf[x) and to perform a modular transformation x — )• 1/x 
in the last integral. One can then choose e = l/vr in order to obtain a numerically rapidly 
converging integral. As a curiosity we observe that Eq.(lO) has a solution for vanishing 
external field i7 = 0, but a non-zero microscopic field B at eB/m^ sa exp(37r/a). This 
would mean that a spontaneous magnetic field would be generated at the Landau-ghost 
pole and therefore lead to a breakdown of Lorentz invariance. Perturbation theory can, 
however, not be naively extrapolated to such very large magnetic fields so any spontaneous 
vacuum magnetization cannot be concluded from the present calculation. 

The magnetic susceptibility, x = OM/dB^ is a measure of the fiuctuations of the mag- 
netization. If only the thermal part of the free energy is retained one would erroneously 
conclude that x — )• [3] and that there would be no fiuctuations in M in the large B 
limit. 
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Figure 1: The vacuum (dotted line) and thermal parts of the magnetization for a neutral 
Fermi gas. Notice that the thermal contribution saturates for large values of the magnetic held. 

3 The Bose Gas 



The formalism used in the previous section apphes also to scalar-QED. The energy spec- 
trum is now given by En = \/k^ + (2n + l)eB + m^ and the one-particle distributions 
are 
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It is rather straightforward to obtain the vacuum correction [7] 
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to the effective action. Proceeding as in the previous section we find the field independent 
thermal part 
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where Zj^iT, fi) is now the partition function for the field independent ideal charged boson 
gas with particle energy E = \/k'^ + m'^ and fsi'-^) = ^(^)/b (^) + ^(~^)/b (~^)- In the 
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Figure 2: The vacuum part (dotted line) and the modulus of the thermal part of the magne- 
tization for a neutral Bose gas. Notice that the thermal contribution actually is negative, and 
that the thermal part of the magnetic susceptibility is changing sign. 



bosonic case we always have |/i| < m. The field dependent part of the effective action can 
then generally, in analogy to the previous section, be written in the form 
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where ??3[2;,(j'] = 1 + 2 X]/^i <?' cos(2/2;). For weak fields, eB <^ m"^, and for T ~^ m the 
last term in Eq.(17) dominates and can be used to calculate the effective temperature- 
dependent coupling as in the previous section. One then uses the expansion YldLi Ko{xl) — > 
TT I {2x) + 1/2 log x as x — )• 0. Even though the coefficient in front of the logarithmic term 
in this case agrees with an asymptotic renormalization group analysis, the term linear 
in T leads to an effective coupling which is a decreasing function of the temperature. 
The difference between bosonic and fermionic QED, in this respect, can be understood 
from a one-loop calculation of the photon polarization tensor. Since we consider a static 
magnetic field we do not obtain any dominant thermal mass of order e^T^. Eor fermions 
the leading high T term is a logarithm which is related to the UV divergence and the 



coefficient is the same as for the standard /3-function. In the bosonic case there is an IR 
divergence at high T from the Bose-Einstein distribution which generates a term hnear 
in T. One can easily verify that the coefficient from the polarization tensor agrees with 
the one obtained from Eq.(17). 

One can, furthermore, verify that the vacuum contribution Ci dominates over the 
thermal parts of the effective action if the magnetic field is sufficiently large. We have 
numerically evaluated Afyac and M^'^ and the result is presented in Fig. 2. We see that 
the vacuum contribution starts to dominate the magnetization of the system when the 
gas changes from a diamagnetic to a paramagnetic behaviour. 

ACKNOWLEDGEMENT 

One of the authors (B.-S. S.) thanks NFR for providing the financial support. It is 
a pleasure to thank the organizers of the 3rd Workshop on Thermal Field Theories, 
1993, for providing a stimulating atmosphere during which parts of the present work were 
initialized. We also thank A. Sjolander for a stimulating discussion. 



References 

[I] p. Flmfors, P. Liljenberg, D. Persson and B.-S. Skagerstam, "Condensation and Mag- 
netization of the Relativistic Bose Gas", preprint Goteborg ITP 94-12. 

[2] V. Ginzburg, ^^High Energy Gamma Ray Astrophysics''^ (North Holland, Amsterdam, 
I99I); S. L. Shapiro and S. A. Teukolsky, ^^Black Holes, White Dwarfs and Neutron 
Stars, The Physics of Compact Objects " (Wiley, New York, 1983); G. Chanmugam, 
"Magnetic Fields of Degenerate Stars", Ann. Rev. Astron. Astrophys. 30 (1992) 143; 
J. R. P. Angel, "Magnetic White Dwarfs", Ann. Rev. Astron. Astrophys. 16 (1978) 
487; J. D. Landstreet, "Magnetic Fields at the Surfaces of Stars^\ The Astron. As- 
trophys. Rev. 4 (1992) 35; R. Narayan, P. Paczyhski and T. Piran, "Gamma-iiaj 
Bursts as the Death Throes of Massive Binary Stars", Ap. J. Lett. 395 (1992) L83. 

[3] D. Miller and P. S. Ray, ^^Thermodynamics of a Relativistic Fermi Gas in a Strong 
Magnetic Field\ Helv. Phys. Acta 57 (1984) 96 and ^^Timiting Particle Density of a 
Relativistic Fermi Gas in a Magnetic FieH\ Helv. Phys. Acta 57 (1984) 704. 



[4] J. Daicic, N. E. Frankel and V. Kowalenko, " Magnetized Pair Bose Gas: Rela- 
tivistic Superconductor"^ Phys. Rev. Lett. 71 (1993) 1779; J. Daicic, N. E. Erankel, 
R. M. Gailis and V. Kowalenko, "Statistical Mechanics of the Magnetized Pair Quan- 
tum Gases", Phys. Rep. 237 (1994) 63. 

[5] P. Elmfors, D. Persson and B.-S. Skagerstam, "The QED Effective Action at Finite 
Temperature and Density\ Phys. Rev. Lett. 71 (1993) 480 and ^^Real-Time Ther- 
mal Propagators and the QED Effective Action for an External Magnetic Eield\ 
Astroparticle Physics (in press). 

[6] W. Dittrich, ^^ Effective Lagrangians at Einite Temperature", Phys. Rev. D19 (1979) 

2385. 

[7] J. Schwinger, ^^On Gauge Invariance and Vacuum Polarization^\ Phys. Rev. 82 (1951) 
664 and ^^ Particles, Sources and Fields^\ VoL 3 (Addison- Wesley Pub. Co., 1988). 

[8] A. Chodos, D. A. Owen and C. M. Sommerfield, ^^Strong Eield Dependence of the 
Eine Structure Constant'', Phys. Lett. 212B (1988) 491. 



